Self-consistent theory of current and voltage noise in multimode ballistic conductors 
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Electron transport in a self-consistent potential along a ballistic two-terminal conductor has been 
investigated. We have derived general formulas which describe the nonlinear current-voltage charac- 
teristics, differential conductance, and low-frequency current and voltage noise assuming an arbitrary 
distribution function and correlation properties of injected electrons. The analytical results have 
been obtained for a wide range of biases: from equilibrium to high values beyond the linear-response 
regime. The particular case of a three-dimensional Fermi-Dirac injection has been analyzed. We 
show that the Coulomb correlations are manifested in the negative excess voltage noise, i.e., the 
voltage fluctuations under high-field transport conditions can be less than in equilibrium. 



I. INTRODUCTION 

Recently, measurements of nonequilibrium noise have 
emerged as a fundamental tool to obtain information 
on the transport properties and, interactions among 
carriers in mesoscopic systems ElotI Shot-noise suppres- 
sion in ballistic conductors caused by Fermi correlations 
has been studied extensively both theoreticallyuo'El and 
experimentallyEi3'ElEj Within the scattering approach, it 
is usually assumed that the ballistic (phase-coherent) 
conductor is attached to reservoirs (terminals or leads) 
with different chemical potentials. In this approach, the 
mean current in a two-terminal conductor is given by 



/= AE / dsT n (e)[f L (s)-f R (e)}, 



(1) 



where q is the electron charge, fh,n{£) the energy distri- 
bution functions at the left (L) and right (R) reservoirs, 
and 7~ n the transmission probabilities (associated with n 
transverse quantum modes (channels). a The correspond- 
ing current-noise PfM^er at zero-frequency has been ob- 
tained in the formeEl (also see Refs. ||| and 11). 



S ' = ^ E I de { r »(e)[/i(l - /i) + M 1 ~ /*)] 

n 

+ % l (s)[l-% l (e)](f L -f R ) 2 }. (2) 

In Eq. (^) the noise is a combination of the thermal emis- 
sion noise of the reservoirs and of the partition noise ap- 
peared due to the current partitioning between the in- 
coming and outgoing states [scattering on tunneling bar- 
rier(s), elastic scatterer(s), point contact(s), etc.]. Al- 
though in some limits the well-known noise terms can be 
identified (associated with thermal noise or shot noise), 
they cannot be separated in general. Out of equilibrium, 
the noise can manifest itself in a different way, depend- 
ing on the conditions. At low biases U, Eq. (g) gives the 
partition shot noise — the excess noise linear in current 
(bias), which does not vanish at zero temperature. It is 
suppressed below the Poisson 2ql value approximately by 
the factoiu J2 n —%)/ J2 n 1~ n - This type of excess 
noise appears whenever there is a partitioning of current 



(T n 0;1). It vanishes for fully ballistic systems for 
which therej-,iS|-fte, partitioning (7~ n =l) (see experimen- 
tal evidenca3'SBI13) . In the absence of partitioning, the 
excess noise is in general no longer linear in the current. 
The inherent randomness in the emission of carriers from 
the reservoirs is at the origin of this type of nonequilib- 
rium noise. Presumably, it is more pronounced for suffi- 
ciently high biases when /jj <C Jl and the transmission 
dominates in only one direction. In this case, noise for- 
mula @i is simplified ted 



[ f L (l-h)de, 



(3) 



where the summation is taken for open channelsS For 
low electron densities the occupation numbers are small, 
ft, <C 1, and Eq. (0) leads to the Schottky formula 
Si = (q 2 /Trh) J2 n * JfLole = 2ql em , where I em is the 
emission current from the reservoir (vacuum-tube-like 
shot noise). In this low-density limit, shot noise is Pois- 
sonian since the transmission of carriers is uncorrelated. 
The factor (1 — /l) in Eq. (||) introduces the Fermi corre- 
lations among carriers when the occupation numbers are 
not small in respect to 1. This leads to the suppressed 
value of the shot noise. Note the difference between the 
partition shot noise mentioned earlier, and the emission 
shot noise given by Eq. (||). The former persists at zero 
temperature, since it reflects the granularity in charge 
transmission manifested by partitioning, while the lat- 
ter vanishes at T — > 0, because its origin is the thermal 
fluctuations of the occupation numbers in the reservoirs. 
Indeed, at equilibrium the sum of two opposite shot-noise 
terms [Eq. (|3|)] gives the Nyquist formula.t3 

It should be stressed that both Eqs. (g) and (||) are 
not complete, since they ignore Coulomb interactions. 
The electrons are charged entities and, while moving 
along the conductor, they affect the electric potential 
giving rise to inhomogeneityEJ The self-consistent cou- 
pling between the nonhomogeneous electron density and 
potential landscape is very important to adequately de- 
scribe the transport an d i aoj se under nonlinear far-from- 
cquilibrium conditions An interesting question is 
how the self-consistency may affect the current and noise 
formulas (0)-((|)? First, the transmission probabilities 
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T n for both current and noise become functionals-of the 
time- averaged self-consistent potential profile ^£3 Sec- 
ond, in the current-noise formulas (||) and (|J), which re- 
flect only the injected current fluctuations SIi n j , the addi- 
tional terms should appear caused by the current fluctua- 
tions blind induced by the fluctuations of the potential Sip. 
Finding the fluctuations 5<p is a complicated problem in 
general, since they are self-consistently linked to the fluc- 
tuations of the occupation numbers along the conductor. 
(The latter may be expressed through the fluctuations of 
the occupation-clumbers at the terminals, since the sys- 
tem is ballistic!!!) . As a result of this self-consistent cou- 
pling, the long-range Coulomb correLatk>»s appear, which 

iBBE2 It is believed, 



may result in the noise suppression! 
however, that such Coulomb correlations need to be 
taken into account for the description of systems in time- 
dependent external fields, or finite-frequency fluctuation 
spectra in stationary fields, while the zero- frequency fluc=. 
tuations in stationary fields are not affected by them.EI 
We show that this is not always true. In the example 
we address in this paper, the additional terms induced 
by the self-consistent field are of the order of the fluctu- 
ations injected from the leads and cannot be neglected 
even in the zero-frequency limit at time-independent bi- 
ases. Moreover, they can almost completely compensate 
the injected fluctuations up to an arbitrarily small value. 
At the same time, the gauge invariance required for the 
charge conservation is fulfilled. We also found that this 
Coulomb suppression of noise is manifested in the nega- 
tive excess voltage noise. Current or voltage fluctuations 
in equilibrium, described by the fluctuation-dissipation 
theorem, usually increase when an external electric field 
is applied. We have an interesting example when an 
interacting (via Coulomb forces) electron system is less 
noisy at far-from-equilibrium conditions than in equilib- 
rium. For noninteracting systems,-such examples have 
been given by Lesovik and Loosen.EJ 

To support our statement we present a theory of cur- 
rent and voltage fluctuations in a ballistic two-terminal 
conductor in a self-consistent field. (Fig. 1). The calcula- 
tion of the self-consistent fluctuating field is, in general, 
a multidimensional problem which includes the electro- 
static environment. For simplicity, we consider rather 
thick samples that allow us to use one-dimensional plane 
geometry for electrostatics. On the other hand, for wide 
conductors the number of transversal modes is large and 
the semiclassical treatment is sufficient. By assuming 
that there is no current partitioning (T n =l for all the 
transmitting modes), we focus mainly on the nonequilib- 
rium noise caused by thermal emission from the reservoirs 
under the action of the long-range Coulomb correlations 
inside the ballistic region, rather than on the partition 
shot noise. It should be noted that the previous theoreti- 
cal studies have been devoted to ballistic conductors with 
a small number of quantum modes (quantum point con- 
tacts) with the Fermi suppression of the partitkwaj-slipt 
noise (Coulomb correlations have been ignored) uuoE3 

The main results of the present investigation are as 



Potential energy 



Electron 
€/ density 




Left 
lead 



Ballistic 
region 



FIG. 1: Energy diagram determining the potential-barrier 
shape for a ballistic two-terminal conductor at equilibrium 
and under applied bias U. Barrier height $(,, as seen from 
the left lead, decreases with bias. The filled area illustrates 
the nonhomogeneous electron-density distribution. 



follows: We have obtained complete analytical expres- 
sions for the steady-state spatial profiles (carrier den- 
sity, self-consistent field), mean current, and differential 
conductance, as well as the current and voltage noise 
powers in ballistic multimode conductors. The analyt- 
ical results have been obtained for a wide range of bi- 
ases: from equilibrium to high values beyond the linear- 
response regime under the self-consistent-field conditions. 
We assume in our derivations arbitrary distribution func- 
tions and correlation properties of injected electrons in 
order to generalize the model to the practically impor- 
tant cases of nanoscale devices with nonequilibrium elec- 
tron injection, like in a resonant-tunneling-diode emit- 
ter, superlattice emitter, hot-electron emitter, etc. (see, 
e.g., Refs. 23l 24,E3,G6h . The particular case of a three- 



dimensional (3D) Fermi-Dirac injection has also been ad- 
dressed. The obtained results clearly demonstrate that 
both the current and voltage noise can be substantially 
reduced owing to the long-range Coulomb interactions. 
This result is very encouraging from the point of view of 
applications. 

The paper is organized as follows. In Sec. II we in- 
troduce the basic equations describing the space-charge- 
limited (SCL) ballistic transport. In Sec. Ill the self- 
consistent steady-state spatial profiles for the electron 
density, electric field and potential are found for an arbi- 
trary injection distribution. The mean current and con- 
ductance are obtained in Sec. IV. Section V describes a 
general formula which relates both current and voltage 
fluctuations with the fluctuations of the occupation num- 
bers in the leads. The current noise power, suppressed by 
interactions, is compared in Sec. VI with the case when 
interactions are disregarded. The Coulomb and Fermi 
noise-suppression factors are discussed in Sees. VII and 
VIII, respectively, whereas the noise temperature is given 
in Sec. IX. The voltage noise power under a fixed-current 
conditions is derived in Sec. X. The implementation of 
the results for a GaAs ballistic conductor is presented in 
Sec. XI. Finally, Sec. XII summarizes the main contribu- 
tions of the paper, whereas in the Appendixes we present 
mathematical details concerning some derivations. 
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II. BASIC EQUATIONS 

We consider a semiconductor ballistic sample attached 
to plane-parallel leads (Fig. |l|) . In a semiclassical frame- 
work, the electron occupation numbers /(x, k, t) inside 
the ballistic conductor are determined by the electron 
flows from the left and right leads. The distribution of 
carriers is nonhomogeneous along the conductor: their 
concentration is higher near the leads and lower in the 
middle of the sample. The inhomogcncity of the space 
charge disturbs the electrostatic potential in such a way 
that the self-consistent built-in field determines the po- 
tential barrier, at which electrons are either reflected or 
transmitted depending on their energy (Fig. [j]). We ne- 
glect tunneling and quantum reflection, i.e., the trans- 
mission probability is 1 if the electron energy is higher 
than the barrier height, and it is in the opposite case. 
Out of equilibrium, the barrier height is different for the 
left and right lead electrons. If for the left electrons the 
barrier height is $h, for the right electrons the barrier 
height is <!>{, + qU (Fig. [I]). This leads to an asymme- 
try in their contribution to the current: as the bias is 
increased, the barrier for the left electrons progressively 
decreases and the current from the left lead enhances, 
whereas the barrier for the right electrons increases and 
the current from the right lead decreases disappearing at 
all at high biases. 

The occupation numbers are described by the Vlasov 
equation (collision less Boltzmann equation with a self- 
consistent field)BEZl 

/ d ^ hk x d ^ ^d(p d \ ~. k ^ q 
\ dt m dx dx Tidk x J ' 

where m is the electron effective mass, k = (k x , kj_), and 
<p(x,t) is the self-consistent electric potential determined 
by the Poisson equation 

d 2 (p q f dk ~ 

dx* = K Jj2^ f{xXt) - (5) 

Here k is the dielectric permittivity and d is the di- 
mension of a momentum space (the spin variable is ne- 
glected). Since carriers move without collisions, the only 
source of noise arises from the random injection of carri- 
ers from the leads. Thus the boundary conditions at the 
left (L) and right (R) leads are: 

/(o,k,t)| fc .>o = /z(k) + a/ £ (M), 

Mk,i,<o = f R (k) + Sf R (k,t), (6) 
<p{l,t) - <p(Q,t) = U(jt), 

where S/l.r are the stochastic forces inside the leads 
with zero average and given correlation properties, and 
U is the applied bias between x=0 and x=l (the po- 
tential drop inside the leads is neglected). As a con- 
sequence of stochastic injection, the occupation numbers 
f(x, k, t) — f(x,k) + 5f(x,k,t) and hence the potential 



ip(x,t) = tp(x) + 6ip(x,t) fluctuate in time around their 
time-averaged values. 

The leads are assumed to be completely absorptive, 
and the transverse electron momenta are conserved. 
Thus, one can make summing up over the transversal 
states (the summation can be replaced by integration due 
to the assumption of a large number of modes) and in- 
troduce for each longitudinal energy e the (fluctuating) 
occupation factor at a cross section x: 

poo 

n(x,e,t)= f(x,e,£±,t)v±de±, (7) 
Jo 

where e—h 2 k x /(2m), ej_= Ti 2 k\ / (2m), and is the den- 
sity of transverse modes = m/2irh 2 for the 3D case). 
The number of occupied transversal modes is N±_ = nA, 
where A is the cross-sectional area. In the semiclassical 
description applied here for a thick conductor, the num- 
ber of the occupied transversal modes is assumed to be 
large, Nj_ > 1. 

III. SELF-CONSISTENT STEADY-STATE 
SPATIAL PROFILES 

It is advantageous to introduce the mean total longi- 
tudinal energy E as a sum of the kinetic energy e and 
the potential energy [—$(2;)]: 

E = e-<f>(x). (8) 

We shall count off the potential energy from the barrier 
top: = q<p(x) — q(p(xb). Therefore, at the barrier 

position x—Xb, we obtain E—0 for electrons having the in- 
jection kinetic energy equal to the barrier height (for both 
leads) . The boundary values for the potential energy are 
expressed through the barrier height qUb and the applied 




FIG. 2: Potential-energy profile for a ballistic space-charge- 
limited conductor. Electrons with energies E > pass over 
the barrier, while those with E < are reflected back to the 
leads. Shadowed regions illustrate the energy distributions of 
the occupation factors at the leads. 
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bias C/ as $i = $(0) = qUb, and = §>(£) = $l + qU states, the electron density at anjz-. section of the con- 
(Fig. |2|). The solution of Eq. (jij) for the stationary case ductor in terms of the potential $Ji§EZl 
(d/dt=0) gives, after integration over the transversal 



JV($)= / [n L {E + $ L ) + n R (E + <S> R )]v{E + <S>)dE + 2 
Jo 



where Ul,r{E) = J °° fL,R{E,e±)i/±ds± are the occupa- 
tion factors at the leads, v{E) = l/[2irfrv(E)] the density 
of states, and v = ^2E/m the velocity. The first inte- 
gral in Eq. (^|) corresponds to the electrons transmitted 
over the barrier (E > 0), while the second integral is 
referred to the reflected carriers (— $ < E < 0). The lat- 
ter term is doubled, since for each energy there are two 
states with opposite momentum: k x and —k x . Finally, 
9{x) with x — x ~ x b is the Heaviside function that dis- 
tinguishes two classes of the reflected carriers: those at 
the left of the barrier (x < 0) originated from the left 



I l9(-x)nL(E + <f> L ) + 6(x)nR(E + $ R )]v(E + $)dE, (9) 

I 

lead, and those at the right (x > 0) coming form the 
right lead. 

The electron density [Eq. (ft)] can now be substituted 
into the Poisson equation d 2 <fr]dx 2 = (q 2 / k)N($) to find 
the self-consistent potential $. The first integration gives 
the electric-field distribution 

m = = -^(x)]JlvW), (10) 

where 



/ ~ — Tfl 



[n L {E + $ £ ) + n R {E + & R )][v(E + $) - v(E)]dE 



2 f [6(-x)n L (E + $ L ) + 0( X )n R (E + ® R )]v(E + <5>)dE 



(11) 



r 



Integrating Eq. @, one obtains the distribution of the 
potential for both regions x < and x > in an implicit 
form 



i\l-x 



-sgn(x) 



Matching the two branches at x = yields 



'I>£ 



d<I> 



w « Jo x/J^m 



VhTm 



(12) 



(13) 



with h- = h(x < 0) and h+ = h(x > 0). Equation (jl^) 
relates three important parameters: the self-consistent 
barrier height Ub, the applied bias U, and the length of 
the conductor I. -Given any two of them, the third one 
can be calculated^! from Eq. (jl|). In Ref. 18 a similar ex- 
pression was obtained for the Maxwell-Boltzmann injec- 
tion distribution. Here we have generalized it to the case 
of an arbitrary injection distribution profiles Ul(E) and 
n R (E) at the leads. It should be noted that in Eq. Jl3| ) 
the dependence on bias enters not only through the up- 
per limits of the integrals, but also through the functions 
h±. 



If we assume that electrons inside the leads obey the 
equilibrium Fermi-Dirac (FD) distributions (the usual as- 
sumption in other works) from Eq. (^) one obtains the 
occupation numbers 



mk R T 
2ttH 2 



ln{l + exp[(e f -e)/k B T}}, (14) 



where Ef is the Fermi energy at the lead and T is the 
temperature. Thus, for the FD case, the steady-state 
spatial profiles of the potential, electric field, and electron 
density are determined through Eqs. (|9])-(|l3|) by making 
use of the distributions 



n L (E + <i> L ) 
n R (E + * a ) 



E/k B T), 

{E + qU)/k B T\, (15) 



where Af = (kpA/4ir) is the number of transversal modes 
in the degenerate zero-temperature limit, £ = Ef/fcgT 
is the reduced Fermi energy, a = (sp — ^ij/ksT is the 
parameter characterizing the position of the Fermi energy 
with respect to the potential barrier, and J~k is the Fermi- 
Dirac integral of index fc.EZl 
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IV. MEAN CURRENT AND CONDUCTANCE 

The mean ballistic current is found as an integral over 
the occupation numbers for the transmitted (E > 0) car- 
riers from both leads: 



/= t4/ [n L (E + $ L )-n R (E + $ R )} dE. (16) 



It is convenient for future analysis to introduce the 
energy-resolved injection currents 



l L,B, 



(E) 



qA_ 

2nh 



n L , R { E )- 



(17) 



By this definition, the mean current is 



the conductance reduces to the value given by the mul- 
tichannel Landauer formula 



(22) 



where N±(^ a L ) = AnL{Q L ) is the number of open modes 
at the barrier energy. Under this small-bias condi- 
tion, the current-voltage characteristics is linear: / ps 
G eq U . Equation ( |2~l] ) may be viewed as the extension of 
the Landauer formula for the conductance to far-from- 
cquilibrium conditions for interacting electrons in a SCL 
ballistic conductor. 

In the opposite limit of high biases Ub -C U < U cr , 
where U cr is the critical bias under which the ha*rier 
vanishes, the asymptotic formula for the current isEZI 



/ = 



[Il(E + $l)- Ir(E + $ r )] dE. (18) j ARjI( 



Having found the barrier height <f>^ from Eq. (|T^), this 
equation determines the current-voltage characteristics. 
Since we have not assumed that the bias must be small, 
this characteristics is nonlinear in bias in a general case. 
For the FD case, Eq. (|l8|) reduces to 



T q . r k B T 



■Tx{a-V% (19) 



where we have denoted the dimcnsionless bias V — 
qU J (ksT). It is seen that under the ballistic SCL conduc- 
tion, the current is determined by the relative positions 
of the Fermi energies at the leads and the barrier top 
through the parameters a and V. This is in contrast to 
the case of diffusive conductors, in which the current is 
determined by the scattering strength. 

The differential conductance G = dl/dU is obtained 
from Eq. (16) as 



n R ($ R ]}. (20) 



The derivative dUb / dll is calculated in the Appendix [A|. 
With its help the formula for the conductance becomes 



n L ($ L )^+n R ($ R )^ 



(21) 



It is seen that the conductance is a sum of two contribu- 
tions corresponding to the left and right leads. Each 
of them is a product of the conductance unit Go = 
q 2 /(2itK) : the number of the transversal modes for the 
injection energy corresponding to the barrier top, and 
some Coulomb interaction factors determined through 
&l,r given in Appendix These factors depends on 
the whole electron system and cannot be separated for 
the left and right lead electrons. 

At small biases close to equilibrium, by assuming iden- 
tical leads (e.g., FD distributions), we obtain Al, r 
A/2. For this case the interaction factors vanish, and 



Child 



f„I L (E + $ L )VEdE 
f™I L (E + $ L )dE 



, (23) 



where the leading-order term is the Child current 



iid 



4 , 
-kA 
9 V m 



2 $ 



3/2 



qP 



(24) 



The main term oc U 3 ^ 2 is independent of the injection, 
while the second-order term oc U contains information 
on the injection occupation numbers. Equation (^3|) for 
the FD case was presented in Ref. ^9|. 

The asymptotic behavior of the conductance at high 
biases is obtained from Eq. (23) 



^ 3 (//child 

= — T 

2 



1 



2 f™I L (E + $ L ) sfEdE 



^ f™I L (E + $ L ) dE 

(25) 

giving the leading-order term Ga ~ (3/2)(/ chi i d /[/) ~ 
VU for an arbitrary injection. For the FD case, Eq. ( |25| ) 
leads to 



^ 3 /child 




qU ^i(a) 



(26) 



V. SELF-CONSISTENT CURRENT AND 
VOLTAGE FLUCTUATIONS. GENERAL 
FORMULAS 

According to the definition of the potential energy in 
Sec. Ill, the potential fluctuations at any point x are 
given by 5<fr x = q5<p(x) — q5<p(xb)- In the nonstationary 
frame fixed to the barrier top, the fluctuations at the 
barrier position are zero, 5§ Xb — 0, whereas at the leads 
they are: (5$ = and 6§i = 5$ R = 8<&l + q5U , 
where 6U is the fluctuation of the applied bias. 

The current fluctuation is obtained by integrating over 
the energy the fluctuation of the occupation factor Sn(E) 



6 



found torn a linearization of Eq. 
values. EH One obtains 



(|J) around the mean 



51 



[ [6I L (E + <!> L )-6I R (E + <f> R )]dE 
Jo 

-Il($l)S$l+Ir@r)8*r. (27) 



where SIl,r{E) are the energy-resolved injection-current 
fluctuations from each lead. In Eq. (p7j), only the low- 
frequency current fluctuations are considered, i.e., the 
frequencies are below the inverse electron transit time 
between the leads and the displacement current is ne- 
glected. The first integral term is standard, and cor- 
responds to the injected current fluctuation SIi n j. The 
last two terms are the fluctuations induced by the self- 
consistent potential fluctuations, that give rise to the 



J 



1l{E) 



1-Ca 



1.3/2 
-E h_ 

9l v(E + $) - v(E) 

o 



3/2 



d<f> 



/■** v(e + $) 

-2C A I 1 Z ' d*, 



1r{E) = 



-1-Ca 



,3/2 



long-range Coulomb correlations.EH To find those terms, 
we need to obtain <5$ L or, equivalently, the self-consistent 
fluctuations of the barrier height in terms of the injected 
fluctuations SIl,r by solving the Poisson equation. This 
has been done in the Appendix |b| The result is the re- 
lation 



/oo 
j L (E) 5I L {E + <5> L )dE 

/oo 
j R (E)5I R (E + $ R )dE, (28) 



where G is the differential conductance [Eq. (21)] and the 
functions jl,r(E) are determined by 



3/2 



3/2 



-$L < ^ < 

, < -E < oo, 

-$r < E < 
, < E < oo. 



(29) 



In Eq. (|29|), one can distinguish the contributions from 
the left-lead (jl) and right-lead (jr) electrons, as well 
as from the reflected (E < 0) and transmitted (E > 0) 
carriers. All the terms related to the barrier fluctuations 
are proportional to the constant 



C A 



2nhA 



(30) 



where A is the constant previously used to determine 
the conductance G and which has been derived in the 
Appendix |A|. Equations (p8|)-(30) are one of the main 
results of our theory. They relate the self-consistent 
current and voltage fluctuations with the noise source — 
spontaneous fluctuations of the occupation numbers in 
the leads. The transfer functions ^l,r, summarizing the 
interaction effects, show the contribution of each energy 
to the total fluctuations. In the absence of interactions, 
j L (E) = 6(E) and j R (E) = -9(E), i.e., the fluctua- 
tions of all energies above the barrier top are equally 
transmitted. The role of the Coulomb interactions is to 
introduce an inhomogeneity in the energy flux of fluctu- 
ations, by suppressing or enhancing occupation- number 
fluctuations at different energies. Note that the Coulomb 
interactions are pronounced only in the presence of trans- 
port. In equilibrium, Ca = 0, and they are not effective. 



In general, both terms jl and 7^ may contribute to 
the noise. However, at high biases, £/& <C U < U cr , one 
can find that only ^l(E > 0) doaainates, the asymptotic 
expression for which is given byE3 



lt(E) = 



VE-va 



1 



VA 



n L (<S>L) Jo 



dn L (E + $ £ ) 
dE 



(31) 



E dE. (32) 



We shall use these formulas later on to analyze the 
asymptotic limits for other important noise quantities. 

To find the total fluctuations SI or SU, one needs to 
define the correlation properties-of the fluctuations at the 
leads. In general, one can writeEH 



(SI k (E)SI k (E')) = K k (E)(Af)S(E - E'), 



(33) 



where, k — L,R and A/ is the frequency bandwidth (we 
assume the low- frequency limit). For the case of the Pois- 
sonian injection from both leads Kl jR (E) oc Jj, ,r(E). 
More generally, for the non-Poissonian injection, under 
the assumption that the leads are in local equilibrium, 
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one can use the formula 



where 



K LiR (E) = 2k B T G A 



dn 



L,R 



dE 



(34) 



where Go is the unit of conductance. This formula follows 
from the Nyquist theorem (see Appendix Q). 

By applying Eq. (K3) for the FD case, we also obtain 



Kl,r{E) 



2G S 



a (-(E/k B T) ' 



(35) 



where Gs = GqM is the Sharvin conductance, and 
Go = q 2 /(2irh) is the unit of conductance. For further 
noise analysis, we have to specify the conditions imposed 
on the external circuit. We shall consider two cases of 
interest: (i) a voltage controlled circuit (zero external 
impedance) for which SU — and one can find the spec- 
tral density of current fluctuations «S7; and (ii) a current 
controlled circuit (infinite external impedance) for which 
SI = and one can find the spectral density of voltage 
fluctuations Sy on the leads. In both cases, we will show 
that Coulomb interactions play a prominent role in the 
noise suppression. 



VI. COULOMB SUPPRESSED CURRENT 
NOISE 



Let us suppose that the potentials at the leads are 
held fixed and do not fluctuate. This corresponds to the 
case when currents are measured using a zero-impedance 
external circuit. Under the condition 5U = 0, Eq. ( pgj ) 
gives the current-noise spectral density 

/DO 
j 2 L (E)K L (E + ^ L )dE 
-Si 

/oo 
lR (E)K R (E + $ R )dE. (36) 

It is important to highlight that the obtained current- 
noise power [Eq. (|36|)1, that includes Coulomb interac- 
tions, has been obtained for a wide range of biases, rang- 
ing from equilibrium to far-from-equilibrium conditions 
beyond the linear-response regime. Therefore, it de- 
scribes both thermal and shot-noise limits. 

One can verify that Eq. ( |36"|) at the high-bias limit 
Ub <S^U < U cr reduces to 



Tl{E)K L (E + Q> L )dE. 



(37) 



Taking into account Eqs. ( pl[ ) and (|34|), one obtains the 
asymptotic expression for the current-noise power 

Sf^/32qI^ = ^ik B TG A , (38) 
qU S 



(3 = 9 



S~I L {E + S> L ) 



dE ' 
2\[E. 



\ 



\ 



Il{<S>l) f™I L (E + $ L )dE 



(39) 



The parameter (3 is determined by the energy profile 
of the injected electrons Il{E). For the EB injection, 
Eq. (p9|) leads to the formula derived earlierE3: 



(3(a) =91- 



7T [Jy 2 (a)] 2 
4^ (a)^i(a) 



(40) 



Equation (g^) is more general and can be applied to an 
arbitrary injection distribution obeying the Nyquist rela- 
tionship [Eq. (34)] for the correlation function. It is seen 
also from Eq. (|38|), that at high biases Sf ~ y/U. 

One can also find, for comparison, the current-noise 
power for the case of disregarded Coulomb correlations 

/•oo />oo 

S™= / K L (E + $ L )dE+ / K R (E + $ R )dE, 
Jo Jo 

which under the assumption of equilibrium conditions at 
the leads [Eq. (|3)] results in 



2qk B T [I l {<Pl) + Ir{®r) 



(41) 



For the sake of completeness, we present also expression 
for the FD case: 

Sj nc °* = 2k B TG s - [Foia) + T Q {a - V)}. (42) 

Note that Eq. corresponds to Eq. (||) discussed in 
Sec. I. Indeed, if one applies Eq. (||) for two opposite 
flows of noninteracting FD electrons, summing up over 
the open channels, one then gets Eq. ((42|). 



VII. COULOMB NOISE-SUPPRESSION 
FACTOR 

To estimate the significance of Coulomb interac- 
tions, |-ane can introduce the Coulomb noise-suppression 
factor^ 



Si 



(43) 



that extends over both thermal-noise and shot-noise lim- 
its. Strictly in equilibrium, Tc = 1, as was pointed out 
in Sec. [V| The effect of interactions is noticeable, how- 
ever, already under small applied biases. In Sec. XI A, 
we will show that while Sj ncor increases with bias, the 
behavior of Si is just the opposite: it decreases with bias 
starting from U =0 up to a certain bias where it reaches 
the noise minimum, then Si increases but much slower 
than Sf DCOT . 
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We remark the difference between the noise-suppres- 
sion factor [Eq. (EJ)] and the s/iot-noise-suppression fac- 
tor (also referred to as the Fano factoru) given by 



F = 



Sj x 
2ql' 



a higher potential barrier by an additional amount of V 
to overcome. It is clear that starting from a certain bias, 
the contribution from only one injecting lead dominates. 
The unidirectional charge flow occurs when 



(44) 



V > max{5; a + 3}. 



(49) 



In the latter formula, the noise power Sf x refers to 
the shot-noise power, pi. e., the excess to the thermal- 
equilibrium-noise levelJa Sj x = Si — Ak B TG eq . Moreover, 
Eq. ( ff~I| ) is meaningful for systems in which Sj x oc I (for 
instance, in linear-response regime). In this case, it sim- 
ply gives a measure of how much the noise power deviates 
from the ideal Poissonian 2ql value due to correlations 
among carriers. For the nonlinear case, when Sf x is not 
proportional to /, definition (El) is less useful, since the 
suppression factor depends on /. It should be noted that 
Eqs. d43| ) and (Q) become identical under the conditions: 
qU S> k B T (for negligible thermal-noise contribution) 
and 5' J ncor = 2ql. The latter is valid, for instance, for 
the Maxwell-Boltzmann nondegenerate injection.Ej If the 
injection is non-Poissonian, as in the case of FD injection, 
S7 ncor ^ 2ql, and Eqs. (EJ) and (El differ. 



This condition comes from the consideration of two lim- 
its. For a nondegenerate injection (the Fermi energy is 
below the barrier, a < —3) the bias-to-temperature ra- 
tio should be large: qU > 5k B T, whereas for a highly 
degenerate injection (the Fermi energy is above the bar- 
rier, a > 3), the bias should be compared with the Fermi 
energy: qU > ef — qUb + 3k B T. Thus, under conditio* 
(|49|) , from Eq. (B8h one obtains the asymptotic formulaEa 



a + 3ir/a 2 ' 



< 



>3. 



(50) 



It is also of interest to analyze the case opposite to con- 
dition (ff9|), when the bias is not so high that both leads 
contribute to the charge flow, namely, V < max{5; a+3}. 
Then, for a nondegenerate limit r F ~ 1- For a highly de- 



VIII. FERMI NOISE-SUPPRESSION FACTOR generate limit, that happens when V < a - 3, one 



It is instructive to introduce the Poissonian noise power 
for the full range of biases: 



use the approximations for the Fcrmi-Dirac intcgralscil 
F Q {x) « x and F x {x) « x 2 /2 + tt 2 /6, and Eq. (EJ) leads 
to a simple formula 



2ql coth 



qU 



2k B T, 
'4k B TG eq , qU^k B T 



(45) 



2 , (V 
r F = -tanh(- 



(51) 



2ql, 



qU » k B T. 



Based on this definition, one can introduce the Fermi 
noise-suppression factor 



T F = 



Sf 



Thus the total noise-suppression factor is 

Si_ 



r = r^rc 



(46) 



(47) 



Note that all definitions Q, (|§), and @ extend from 
thermal- to shot-noise limits. They will be used in the 
analysis of the results in Sec. XI. 

For the FD two-lead injection, from Eqs. ( |l9| ) and (EJ) 
it follows that 



Surprisingly, the dependence on a and hence on the bar- 
rier height, ballistic length, and material parameters, 
canceled out from this equation. Tp depends only on one 
parameter — the bias-to-temperature ratio, and at suffi- 
ciently high V it decreases with bias as Tp ps 2/V. Note 
that this behavior occurs under the nonlinear bias regime 
in the presence of a space charge. Indeed, for this case one 
finds the sublinear characteristics for the mean current 
I = GsU[l — q(Ub + U/2)/ep], whereas the current noise 
power is given by Sf nml = 4k B TG s [l - q(U b + U/2)/e F ), 
with the identical factor in square brackets. As a re- 
sult, one obtains, Sj ncor — 4k B T I /U, as in the linear- 
response regime. It should be remembered, however, that 
at higher biases V ^S> a, the 2/V dependence is changed 
to the 2/ a law. The largest noise suppression by Fermi 
correlationscj is described by Eq. ( |50| ) giving Tp m = 2/£, 
when the barrier height is zero. We shall give some ex- 
amples in Sec. XI. 



f — 



tanh(V/2). 



(48) 



One can verify that at the low-bias limit V — * 0, there 
is no suppression effect: Tp — > 1. The finite bias intro- 
duces asymmetry in the contributions from electrons of 
different leads. The larger the bias, the smaller is the con- 
tribution from the biased lead, since the electrons have 



IX. NOISE TEMPERATURE 

It is interesting to see from Eq. (|38|), that at high bi- 
ases and strong screening, despite the strong nonlinearity, 
the ratio between Si and the differential conductance G 
tends to the constant value. It is instructive, therefore, 
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to introduce the effective noise temperature T n throughc2l 
Sj = A)tBT n G. Note that ksT n has a meaning of the 
maximum noise power per unit bandwidth which can 
be delivered to an output matched circuit, thus it is a 
measurable quantity.Eil The asymptotic high-bias value 
is then obtained as 



T 



1 



/3(a) 



3(l-7r/4), a 



< 



1/3, 



a > 1. 



(52) 



It is seen that T n < T for any a, indicating the noise 
suppression effect [see the plot of /3(a) in Ref. ^9|. For 
a nondegenerate Maxwell-Boltzmann injection, the lim- 
iting vakie, (T^ d °s/T) = 3(1 - tt/4) ps 0.644 is well 
known .Mr 2 ! For a highly degenerate FD injection, we 
have obtained from our theory (T r f s /T) = 1/3. The 
physical meaning of the latter result is that the noise 
power per unit bandwidth produced by the SCL ballis- 
tic conductor with degenerate FD electrons is 1/3 of the 
thermal noise power produced by the heated resistance 
with the same value of the conductance G (the same I-V 
dependence), independently of the material parameters. 



X. COULOMB SUPPRESSED VOLTAGE NOISE 

Alternatively, one could measure the voltages at the 
leads using an ideal infinite-impedance voltmeter. The 
infinite-impedance external circuit then forces the cur- 
rent to be zero at all times, 81 = 0. Fluctuations in the 
current are counterbalanced by fluctuations of the chem- 
ical potentials in the electron reservoirs. Under the con- 
dition 61 = 0, Eq. ( |28| ) gives the voltage- noise spectral 
density which takes into account the Coulomb correla- 
tions: 



Su 



1 

G 2 



/OO 
7 2 L (E)K L (E + <f> L )dE 
-*z, 



lR {E)K R (E + $ R )dE 



It is evident that the relation 



Si 



5j 
G 2 



(53) 



(54) 



holds, in which Si is the current-noise spectral density 
[Eq. (p6^] measured under 6U = 0, and G is the steady- 
state differential conductance (pl|). 

The asymptotic behavior of Su at high biases can also 
be found from Eqs. ( |38| ) and (|5J). We obtain 



S£ = ^4k R T- 



G A 



(55) 



with /3 given by Eq. (|39|). It is seen that the voltage noise 
decreases with bias as S$ ~ 1/ y/U at t/& <C U < U cr (the 
general result independent of the injection distribution). 
Hence the Coulomb interactions result in the voltage- 
noise suppression. This fact will be discussed in detail in 



Sec. XI B. 



XI. RESULTS FOR FERMI-DIRAC INJECTION 

To illustrate our results— consider the GaAs ballis- 
tic n-i-n diode at T=4 Ko For this temperature and 
the effective mass m=0.067mo, the effective density of 
states is N c rj 6.7 x 10 14 cm~ 3 . Assuming the con- 
tact doping 1.6 x 10 16 cm -3 , the reduced Fermi en- 
ergy £ ps 10, and the contact electrons are degener- 
ate, that is necessary for studying the joint effect of 
both Fermi and Coulomb correlations. For this set of 
parameters, the Debye screening length associated with 
the contact degenerate electron density is approximately 
L D = ^/Kk B T/{q 2 N c F_ 1/2 (Z)} fa 14 nm. The calcu- 
lations have been carried out for the following ballistic 
lengths: £=0.05, 0.1, and 0.5 ^m. 

The degeneracy of the contact electrons does not guar- 
antee the degeneracy of the injection, since the poten- 
tial barrier determines the energy portion of electrons 
which may pass over the barrier and contribute to the 
injection current. For each ballistic length £ and bias 
[7, we have solved numerically Eq. (|l3|) to find the self- 
consistent potential barrier height <S>l and the parameter 
a = (ep — $i)/fcsT characterizing the position of the 
Fermi energy with respect to the potential barrier. The 
results are plotted in Fig. [| 

It is seen that in equilibrium, for the ballistic lengths 
£=0.05, 0.1, and 0.5 /im, the self-consistent barrier height 
<f>L is about 6 fc^T, 9 fcsT, and 13 k R T, respectively. 
This means that the injected electrons at U=Q are de- 
generate for £=0.05 and 0.1 /im, and nondegenerate for 
£=0.5 /im, since for the latter case only the tail of the 
distribution function is injected (a < —3). As U is in- 
creased, $l vanishes and the injection becomes degener- 
ate for all three cases. Finally, at U — > U cr the poten- 
tial barrier vanishes, a=£, and the transport is no longer 
space-charge limited. The values of U cr depend obviously 
on the ballistic length £ (see Fig. ||). 



A. Fixed-bias conditions: SV = 

In Fig. U we show the results for the current-noise 
power Si [Eq. (|36j)] and differential conductance G vs 
bias U calculated for different ballistic lengths I. For 
comparison, the noise power for the Poissonian injection 
Sf [Eq. (|45|)] and Fermi-Dirac injection with disregarded 
Coulomb correlations Sp 11 ™ 1 ' [Eq. (|4^)] have also been 
plotted. In equilibrium, all the noise-power curves coa- 
lesce toward the Johnson-Nyquist noise Sj q = AksTGeq 
independently of the presence of Fermi or Coulomb cor- 
relations. However, starting from small biases the differ- 
ence becomes drastic. While Sf ncor increases with bias, 
the behavior of Si is just the opposite: it decreases with 
bias starting from U—0 up to a certain bias where it 
reaches the noise minimum, then Si increases, but much 
slower than Sf ncor . Finally, at U — > U cr when the barrier 
vanishes, Si sharply recovers 5} lncor . Note that in the ab- 
sence of Coulomb interactions, Sj ncor follows the Poisso- 
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FIG. 3: Parameter a = (ef — $i)/fcsT characterizing the 
position of the Fermi energy ef with respect to the poten- 
tial barrier <&l for different lengths of the ballistic region I. 
At U — > U cr the barrier vanishes, $l=0, and a attains its 
maximum value a=£=10. 



T n < T would indicate the significance of the Coulomb 
correlations effects which suppress the current noise. At 
U ~ U cr , T n sharply increases due to the current satura- 
tion (G=0), that may also be detected in the experiment. 

The current-noise-suppression factors Fq and Tp, and 
the total r (their multiplication) correspondent to the 
noise-power curves of Fig. |i| are plotted in Fig. [?] as 
functions of bias. The behavior of the Fermi suppres- 
sion factor Tp is in a clo se agreement with the analyti- 
cal formulas of Sec. [VIH . It varies from 1 at low biases, 
then decreases attaining with a nice precision the asymp- 
totic dependence 2/(a + 3n/a 2 ) at high biases [Eq. (ft7|)]. 
For all three cases, the same minimal value 2/£ w 0.2 
is reached at, the highest biases, in agreement with the 
predictions .Ea We have also checked that for the degen- 
erate injection from both leads, which is well realized for 
the length £=0.05 /zm, the analytical formula [Eq. Q)] 
very nicely describes the numerical results in a wide bias 
range : from U=0 up to U ~ 5k B T [Fig. @(a)]. 

The relative significance of two mechanisms on the 



nian law Sf = 2ql coth(qU/2ksT) only for a nondegen- 
erate injection, as in the case of t =0.5 /zm at qU < lOfc^T 
[Fig. f|(c)]. At higher biases, and for shorter ballistic 
lengths in all the range, the injection is non-Maxwellian 
and Si < Sf because of Fermi suppression [Figs. H(a) 
and §(b)] . It should also be noted that Si < Sj in 
a wide range of biases, which means that the noise for 
interacting ballistic electrons in an external field is less 
than the equilibrium Johnson-Nyquist noise. As will be 
shown later, the same is true for the voltage noise power. 

Figure || assists the understanding of the results by 
showing the contributions to Si from different electron 
groups. The effect of Coulomb correlations is manifested 
quite differently for the left-lead and right-lead electrons: 
while the left-lead noise is suppressed (SiLt < <5y£ cor ), 
the right-lead noise is enhanced (Sim > Sf^^). Since 
the role of the right- lead electrons is diminished with bias, 
the overall effect of interaction results in the total-noise 
suppression. There is also a non-negligible contribution 
(~ 10 — 15%) to the noise from the reflected carriers at 
qU < lOfcsT. It appears for correlated electrons only. 

Figure ^ shows the noise temperature T n versus bias 
U calculated from (T n /T) = Si/(4k B TG) by using the 
data of Fig. [|. One can see that starting from T n =T at 
zero bias it drops at qU > k B T below the temperature 
T of the injected electrons. It is interesting to note that 
for degenerate electrons (£=0.05 /im), this drop starts to 
appear at higher biases than for nondegenerate electrons 
(£=0.5 ^m). According to Eq. (|52|), the minimal asymp- 
totic value of T n in the limit £ — > oo, U — > oo differs 
for nondegenerate and degenerate electrons (see indica- 
tions in Fig. ||). For our set of parameters, the injec- 
tion is degenerate at the highest biases. However, the 
limit T T f 9 = (1/3)T is not achieved for those ballistic 
lengths, since the samples are not sufficiently long. Note 
that T n is a measurable quantity, and the observation of 
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FIG. 4: Current-noise power Si and differential conduc- 
tance G vs bias U for different ballistic lengths I. For 
comparison the noise power for the Poissonian injection 
2ql coth(gi7/2fcgT) and the Fermi-Dirac injection with dis- 
regarded Coulomb correlations gy ncor are plotted. The nor- 
malization constants are the corresponding equilibrium values 
at [7=0: Sj q = 4k B TG eq and G eq . 
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noise suppression can be understood by comparing the 
curves for three different ballistic lengths £. For short 
samples (£=0.05 /im), the Fermi suppression dominates 
at low biases V < 3, where ~ 1 and T w [see 
Fig. 0(a)]. At higher biases, 3 < V < 20, both Coulomb 
and Fermi mechanisms contribute to the suppression. For 
£=0.1 jiim, the suppression factors Tp and Tc are com- 
parable in all the bias range. Finally, for the longer sam- 
ple, £=0.5 /urn, the Coulomb noise suppression completely 
dominates: Tc <C Tp. This behavior can be explained 
by the fact that the Fermi shot-noise-suppression factor 
is limited below by the value 2/£, i.e., by the proper- 
ties of the injecting contact independently of the ballistic 
lcngth.Ea In contrast, the Coulomb noise suppression may 
be enhanced arbitrarily strong by extending the length 
of the ballistic sample with a simultaneous increase of 
bias (provided the transport remains ballistic). There- 
fore, for any degree of electron degeneracy £ = ef/AsbT 1 , 
there exists the ballistic length starting from which the 
Coulomb interactions become to dominate in the noise 
suppression. 

Another important difference between the two suppres- 
sion mechanisms is that the Fermi noise-suppression fac- 
tor IV is a monotonically decreasing function of bias, 
while the Coulomb noise-suppression factor Tc exhibits 
a minimum at a certain bias value, as seen in Fig. |t]. Af- 
ter the minimum, the curve of Tc increases to 1 due to 
the disappearance of the potential barrier at U = U cr . 

The total noise-suppression factor T approaches at 
high biases the asymptotic curve fc^T '/qU , once the in- 
jection, because of barrier lowering, becomes fully degen- 
erate, in agreement with the prediction.Ej The longer the 
sample, the wider is the bias range in which this asymp- 
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FIG. 5: Contributions to the current-noise power Si for the 
case of £=0.5/im: from overbarrier electrons transmitted from 
the left (Si u) and right (Sim) leads, and those reflected 
by the barrier (Si Lr+Rr)- For comparison, contributions to 
ST cor are shown as well: from left- (5j| cor ) and right-lead 
(Sr R 01 ) electrons. 
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FIG. 6: Noise temperature T n vs bias U for different bal- 
listic lengths £. The asymptotic (I — > oo, U — > oo) limits 
for nondegenerate and highly degenerate electrons are shown: 
T" dea = 3(1 - tt/4)T and T^ 9 = (1/3)T, respectively. 
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FIG. 7: Current-noise-suppression factors Tc (Coulomb), Tf 
(Fermi), and r=IVrc (total) as functions of applied bias U 
for different ballistic lengths I. The analytical approximations 
2/(q + tt 2 /3q) and 2/Ftanh(F/2) for I> and k B T/qU for T 
are also shown. 
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FIG. 8: Energy distributions of the current-noise power 
s(E) under various biases V for a ballistic conductor with 
fc0.05^tm (solid lines). The zero energy corresponds to the 
conduction-band edge at the left (unbiased) lead. The sharp 
peak at low energies corresponds to the position of the space- 
charge barrier. The dashed line shows the energy profile 
Kl(E) for the injection noise at the left lead (ef = lOfcsT). 
The profile for the right lead Kr(E) is the same, but is shifted 
by —V in energy. All the curves are normalized by iqlL^I) 
related to the noise level at V=0. 



totic law is fulfilled independently of the material param- 
eters (Fig. 0). It is also important that the suppression 
may be several orders of magnitude stronger than the 
shot-noise suppression due to elastic partitioning.tj 

It is instructive to plot the energy-resolved current- 
noise power s(E) defined by Si — J s(E)dE. The de- 
rived formula [Eq. ([36])] allows us to analyze these dis- 
tributions for different lengths and biases. The results 
for <?=0.05 and are shown in Figs. || and |[ At 

small biases (V^=0.1 in the figures), the Coulomb interac- 
tions are ineffective, and the noise is approximately the 
sum of two equal contributions from the left and right 
leads. These contributions are the Fermi-Dirac profiles 
filled out above the barrier (the contributions of the re- 
flected carriers with energies below the barrier is negli- 
gible at these biases). With increasing the bias, several 
features appear: (i) the contribution from the right con- 
tact becomes smaller and smaller because of the shift in 
energy — V, (ii) the Coulomb interactions give rise to a 
sharp peak at the barrier energy with a noise suppression 
at the energies beyond the peak, and (iii) the carriers be- 
low the barrier give appreciable nonzero contribution to 
the noise. The peak appears due to the fact that electrons 
with the energy E = virtually stop at the barrier top, 
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FIG. 9: Distributions similar to Fig. ^ for another ballistic 
length ^=0.5^im. The notations are the same, except that 
Kl(E) is not normalized. 

producing an infinitely large perturbation of the current 
(this singularity is integrable, since it is of the logarith- 
mic typelia). Another interesting feature is the "noise- 
less" energy E* lying above the barrier, in which the 
noise exhibits a local minimum. It is better pronounced 
for a nondegenerate injection (see, for instance, Fig. ||) 
where one can observe the zero- noise point s(E*)=0 for 
various biases) . This point appears approximately at the 
condition: jl(E*)=0. As long as the barrier vanishes at 
highest biases U — > U cr , the Coulomb noise suppression 
disappears and the energy profile s(E) recovers the FD 
shape. 

B. Fixed-current conditions: SI = 

Thus far we have presented the results obtained under 
the assumption that the ballistic sample is connected to 
zero-impedance external circuit. In this case the fluctu- 
ations of the applied voltage can be neglected. In exper- 
iments, it is the voltage fluctuations which are actually 
measured and which eventually are converted to current 
fluctuations. By using an infinite-impedance circuit, the 
current fluctuations are forced to be zero, and one can 
analyze the voltage-noise power. Both cases are interre- 
lated through Eq. (|4]). It is of interest however to see 
the results for the voltage-noise power - the quantity that 
can be measured directly. 

Figure [To] shows the results of applying of Eq. ( |54"| ) to 
our set of parameters. The behavior of Su calculated 
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with and without Coulomb correlations is strikingly dif- 
ferent. We remark the following features: (i) For the case 
when the interactions are included, the noise decreases 
with bias instead of increasing [the asymptotic behavior 
Sjj ~ 1/VU at high biases (see Sec. |x|) is confirmed], 
(ii) For longer samples, the range of the space-charge 
conduction is wider, and the suppression of voltage fluc- 
tuations is much more pronounced, (iii) Comparing the 
asymptotic dependences Sjj ~ 1/ VU and Sj ~ y/U, it 
is seen that the latter eventually exceeds the equilibrium 
Nyquist noise when the ballistic length is sufficiently long 
[see Fig. ||(c)]. In contrast, Su falls off the equilibrium 
value in a full range of SCL conduction. 
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FIG. 10: Voltage-noise powers Su vs bias U for different bal- 
listic lengths £ are compared for two cases: with and without 
Coulomb correlations taken into account. 

Note that the equilibrium fluctuations, described 
by the fluctuation-dissipation theorem, usually increase 
when an external electric field is applied. In our model, 
we have an interesting example when an interacting (via 
Coulomb forces) electron system is less noisy at far-from- 
cquilibrium conditions than in equilibrium. For noninter- 
acting electrons in quantum conductors, such examples 
have been provided by Lesovik and Loosen.EJ 

We can also mention more familiar examples from 
semiconductor device literature. For instance, in 
Schottky-barrier diodes or p-n junctions, in the range 
of the exponential I-V characteristics, the current-noise 
power is given by Sj = 2IzbTG, which is a half of the. 
thermal noise value given by the Nyquist relationship.til 
However, this is not really a suppression effect, since the 
current-noise power Si never drops down the equilibrium 
Nyquist level Sj q — 4kaTG e q- In this case G 3> G eq and 
Si > Sj 9 for any biasE3 (in our model Si < S^ q in a 
wide bias range). On the other hand, for these junctions 
under the fixed-current conditions, Su = 2k,BTG~ 1 , that 
is again a half of-tlie Nyquist relationship, but Sjj < S v q 
may now occur It should be emphasized, however, 
that the latter noise reduction appears for noninteract- 
ing carriers, and it is caused by the nonlinearity in the 
current- voltage characteristics which results in such a be- 
havior that the conductance G grows with bias as fast as 



the current-noise power Si (exponentially). As a result, 
Su = Si/G 2 is a decreasing function of bias. In our 
model, the noise suppression below the thermal equilib- 
rium value (negative excess voltage noise) occurs due to 
Coulomb interactions among carriers. Without interac- 
tions, despite the nonlinear SCL regime, the noise grows 
above the Nyquist level, as was shown in Figs. |] and [!(]. 



XII. DISCUSSION 

In this paper, we have presented a theory of the elec- 
tron transport and noise in a self-consistent potential 
along a ballistic two-terminal conductor. Since electrons 
are fermions and carry charge, they interact among them- 
selves by both Fermi statistical correlations and long- 
range Coulomb correlations. The interplay of these two 
mechanisms determine the noise properties of a ballistic 
conductor, - the subject we have addressed in the paper. 

The long-range Coulomb correlations appear due to 
the self-consistent coupling between the electric poten- 
tial and the occupation numbers. This coupling is essen- 
tial to adequately describe the noise phenomena. To de- 
velop a better understanding of the Coulomb-correlation 
effect, we rewrite Eq. (|l|) — a standard equation for the 
mean current in a two-terminal conductor — in which we 
explicitly introduce the dependence of the transmission 
probabilities on the self-consistent potential (in this case 
on the barrier height <&&): 



1=1 de[T L (E,$ b )I L (s)-T R (s,$ b )I R (£ 



)]. (56) 



Here, Ilm{ £ ) = {qA/2-Kn) J fiM £ > £±)v±de± are the 
currents corresponding to the longitudinal energy e. 
In the semiclassical limit, by neglecting the quantum- 
mechanical reflection of electrons with energies e > $b 
and tunneling through the barrier, the transmission prob- 
abilities are the Heaviside step functions: Tu (e,$b) = 
9{e - $ b ) and T R (e, $ 6 ) = 9{e - $ b - qU). 

The current fluctuation is found by perturbing 
Eq. (K9); we obtain 



51 = J[T L SI L -T R SI R ] ale 



= SL, 



(S* b ) (S<f b - qSU) ^ I R 

oe oe 

•A "T" SIi n d^ 



de 
(57) 



The first integral in Eq. (57) is the injected current fluc- 
tuation SIi n j . Its origin is the thermal fluctuation of the 
occupation numbers in the leads ballistically injected into 
the conductor. In fact this term corresponds to Eq. (||) - 
the standard formula used to calculate the current noise 
in mesoscopic conductors. [Since for our case all T n = 
or 1, the term oc T n (l — T n ) is absent.] The second inte- 
gral in Eq. ([57]) is the induced current fluctuation SLmd 
caused by the fluctuation of the potential. It is precisely 
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the term appeared due to Coulomb correlations and ig- 
nor ed in Eq. (§). For our case, the derivatives are found 
as: (dT L /de) = S(e - $ h ), (dT R /de) = 5(e - * 6 - qU), 
leading to 

SI md = -S* b I L ($ b ) - (S$b - qSU) I R ($ b + qU). (58) 



Thus Eqs. (J57|) and (|58|) lead to Eq. (|27j) for the cur- 
rent fluctuation derived more rigorously earlier from the 
transport equation. 

We would like to highlight that the induced current 
fluctuations 5Ii n d should appear not only in the case of 
completely open/closed channels (T n — 0;1), but also 
under the conditions of the partitioning shot noise, for 
which there exist channels with < T n < 1. It is clear 
that Blind should depend in general on the derivatives 
of the transmission probabilities (dT n /de) and the fluc- 
tuations of the self-consistent potential 5$> x . The main 
problem is then to find the fluctuations 6& x through the 
noise sources. For the particular case of a multimode 
ballistic two-terminal conductor, we have found an exact 
analytical result for SQ X . For the case of partitioning 
shot noise in which Eq. (|J) holds, work is in progress. 

The validity of our theory can be tested experimentally 
in currently accessible semicoftductor structures. The re- 
quired conditions are similartEl to those for the transport 
in vacuum tubes: (i) the ballistic electron transmission 
between the terminals, and (ii) the limitation of current 
by the space charge. The SCL transport regime, as ap- 
plied for ballistic electrons in solids (mostly in n + -n-n + 
or n-i-n semiconductor strictures), was discussed a long 
time ago (see, e.g., theoryE3 and experiments^) . Unfor- 
tunately, the data, on noise measurements in these struc- 
tures are scarce. Ej Due to a great progress in noise mea- 
surements. m,ouantum ballistic conductors during the last 
ten yearsLMmLO (also see Ref. ^), we believe it would be 
now possible to measure the noise suppression effects in 
SCL ballistic conductors. Although the theoretical re- 
sults presented in this paper are strictly valid for thick 
multichannel conductors (3D electron gas) , the Coulomb 
suppression of noise should also be pronouncedB in con- 
ductors with a small number of channels (2D or ID) in 
which electrons are more confined in space, for instance, 
in quantum Jttircs under the high-bias nonlinear trans- 
port regimej^l or in carbon nanotubes under the SCL 
conduction!^ 

Additionally, we would like to emphasize the impor- 
tance of the effect of Coulomb interactions. They not 
only lead to the noise reduction, but can also be used as 
a tool to probe the energy profile, of the injected carriers 
and other electronic properties.CJ 
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APPENDIX A: DERIVATION OF dU b /dU 
Differentiating Eq. (|l3|) gives 



dt>, 



du VmH) 2 j h m du 



1 dh_ 



+ 



d$. 



dU 7M$^) 2 J Q ,3/2 dU 



1 dh 4 



d$ = 0. (Al) 



By using d$ R /dU = q+d^L/dU and finding dh/dll from 
Eq. (pdl), we obtain 



dh 
dU 



-(H L + H R )^-qH R , (A2) 



where 



m 

2^h 



jjf T> L [v{E + 3>)-v(E)}dE 
+29(- x ) J V L v(E + *) dE j , (A3) 
H R {$) = — |y V R [v(E + $) - v(E)] dE 

+29( X ) J _V R v(E + $) dE J , (A4) 
, k = L,R. (A5) 



V k (E) 



dn k 
dE 



Substituting Eqs. (A2)-(A4) into (Al), we finally obtain 



dU b ld<f> L 



dU q dU A L + A R 
Here we have defined 



(A6) 



A R = 



y/h-($ L ) ' Jo h 3 / 2 



?R m 



o h 



3/2 



o hi 



3/2 



where we have denoted HZ = H/-(x < 0) and H£ 
Hi ; \ > 0), k = L. R. The quantities Al and A R , as 
well as their sum 



A = At, + A f 



(A7) 



are used in this paper to calculate the differential conduc- 
tance [Eq. (f2l|)] and the noise suppression [Eq. (p3C|)]. The 
physical meaning of A becomes clear from the relation 



A n n ( 3i 

a = 2 "V«fe [ - 



(AS) 



i.e., it relates the increment of the barrier height with 
the increase of the length of the sample under fixed bias. 
Ah.R are the corresponding contributions to that incre- 
ment from the left-lead and right-lead electrons. 
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APPENDIX B: DERIVATION OF THE 
SELF-CONSISTENT POTENTIAL 
FLUCTUATIONS 



electron-density fluctuation as a sum of two contribu- 
tions, SN = 8N ln: > + SN md , where the injected part is 



Integrating the fluctuation of the occupation factor 
8n(E) over the longitudinal states, one obtains the 



J 



5N 



rOO 

'"■?($) = / [6n L (E + $ L ) + 6n R (E + $ R )]v{E + $)dE 
Jo 

+ 2 f [0(-x)6n L {E + $L) + 6(x)SnR(E + $ R )]v(E + $)dE, (Bl) 



r 



and the induced part is 

SNind W = ~ ^T»* ^ST«*- ( B2 ) 

c?<f> d$ d$ 



Equations (Bl) and (B2) should now be substituted into 
the Poisson equation for S$ x to find the self-consistent 
fluctuation of the potential profile: 



d 2 a 2 

S<^ x = — (5N inj + 5N md ). 
dx 2 K 



We get 



dx 2 



q 2 dN 
n d<& 



5$ x = Ss x , 



(B3) 



(B4) 



where Ss x 
(dH R /d$)5<f> R 



(q 2 /K)[5N ln i - (dH L /d<Z>)S<f> L - 
The boundary conditions for this 
equation 6® x (0) = <5$ L , = 6$> R , 5<£> x (x b ) = 0. 

The second-order differential equation (B4) with spa- 
tially—dependent coefficients can be solved explicitly for 
<5$ X .E£| Here we need just the boundary values S$l and 
6$>r (the relation between them), which has entered ex- 
plicitly into the nonhomogeneous part and can be ob- 
tained by applying the Green's identity for the operator 
L 



[u(x)L5$ x — 8<& x Lu{x)]dx 



(B5) 



where [a; b] = [0;xb] for x < and [a; b] — [xi,;l] for 
X > 0. It is convenient to chose the function u(x) as a 
solution of the homogeneous equation Lu(x)=0 satisfying 
the boundary conditions u(0)=0 and u(l)=0. This gives 



/ u 8s x dx + u 5s x dx 
'o Jx b 



6<f>, 



5<I>, 



where £l and £ R are the electric fields at a;=0 and 
x—£, respectively. Changing the variable of integration 
dx——d^/(q£) , one gets 



£ 



■ 8s x d$ 



£ 



5s x d$ 



5<S> L 6$} 



£l £r 



(B6) 



It is convenient to represent the fluctuation 5s x as a 
derivative 5s x — (n/q 2 ) (d5h/d&). By using this, nota- 
tion, the integrals in Eq. (jB6|) can be reduced tcE3 



k u ddh 



1 



8h 



k = L,R (B7) 



, £ d$ " qJo £ 3 
whereas Sh is obtained by integration of Ss x : 

Sh = Sh" 1 ? - H L 6$ L - H R SQr, (B8) 

5h in3 {<5>) = / 6N inj d$. (B9) 
Jo 



Now substituting Eqs. (gj and into Ec l- (@)> and 
by using Eq. (|l0|), we obtain 



AS^ L +A R qSU = I -jjjd® 



** 5hT 



,3/2 



(BIO) 



where A and Ar were d enoted in Appendix [X]. Combin- 
ing Eqs. (^7|) and ( B10| ) and excluding <5$l, we obtain 



/•oo 

SI-G5U= / [5I L (E + Q L ) -5I R {E + <S> R )}dE 
Jo 

- Ir{$r) ( f* L Sh m > „ /"** ShT 



,3/2 



3/2 



that leads to Eq. 
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APPENDIX C: NYQUIST THEOREM AND THE 
BOUNDARY CONDITIONS FOR 
FLUCTUATIONS 



Since this integral relation should be valid for differ- 
ent lengths I of the ballistic conductor (different $°), 
it should also be valid for the integrands, 



Consider the situation when the potentials at the 
leads are held equal (t/=0, SU=0) by means of a zero- 
impedance external circuit. (A similar consideration can 
be carried out for the infinite-impedance circuit.) Ad- 
ditionally we assume that the contacts are identical: 
Il{E)_= Ir(E), and K L {E) = K R (E), VE. Thus from 
Eq. © we have C* A =0, j L (E) = 6(E), and j R (E) = 
—9(E), which means that Coulomb correlations do not 
affect noise at zero bias. Therefore, from Eqs. (28) and 
(p3|) one obtains the current-noise power 



/>oo 

S? = 2 / K L (E + <f>° L )dE, 
Jo 



(CI) 



where <f>° i s the equilibrium barrier height (the noise de- 
pends on the steady-state self-consistent field). For the 
equilibrium conductance we find 



eq ~ dU 



= -q 

u=o Jo 



dE 



E+<S>°, 



dE. (C2) 



By using the Nyquist theorem Sj q — Ak B TG eq , we obtain 



/ K L {E + ^° L )dE = 2qk B T 
Jo Jo 



d_h 

dE 



dE. 



K L (E) = 2qk B T 



dE 



(C3) 



that leads to Eq. (p4[). Thus, just from the Nyquist theo- 
rem we have a useful relation for the energy-resolved cur- 
rents (occupation factors) at the leads. It relates the en- 
ergy profiles of the fluctuations and the mean values. In 
the simplest case of the Poissonian injection, for instance, 
the correlation function is proportional to the meanO 



K%™(E) = 2qI L (E). 



(C4) 



From this result it follows that h(E) = -k B T(dI L /dE), 
and one obtains the Boltzmann distribution 



h L ™{E) =Cexp(-E/k B T), 



(C5) 



where the integration constant C is determined by the 
normalization condition. 
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